home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 40 / Amiga Format CD40 (1999-05-11)(Future Publishing)(GB)(Track 1 of 3)[!][issue 1999-06].iso / -readerstuff- / paul_qureshi / source / nbody.lzx / NBody_Collision_Detection / project / cylinder.c < prev    next >
C/C++ Source or Header  |  1999-03-27  |  1KB  |  50 lines
  1. #include <stdio.h>
  2. #include <math.h>
  3.  
  4. #define PI 3.141592653589793
  5.  
  6. main()
  7. {
  8.  
  9.     int n, m, i=0;   /* loop counter                  */
  10.     double x, y, z;     /* coordinates              */
  11.     FILE *fpr;
  12.  
  13.     fpr = fopen("cylinder.dat", "w");
  14.     fprintf(fpr, "polyhedron   cylinder\n");
  15.     fprintf(fpr, "n0  0  0  2.01\n");
  16.  
  17.     for(m=0; m<32; m++) {
  18.            x = 2*cos(m*PI/16.0);
  19.            y = 2*sin(m*PI/16.0);
  20.            z = 2;
  21.            fprintf(fpr, "n01.m%02d   %2.14f  %2.14f  %2.14f\n", m,x,y,z);
  22.          }
  23.  
  24.     for(m=0; m<32; m++) {
  25.            x = 2*cos(m*PI/16.0);
  26.            y = 2*sin(m*PI/16.0);
  27.            z = -2;
  28.            fprintf(fpr, "n02.m%02d   %2.14f  %2.14f  %2.14f\n", m,x,y,z);
  29.          }
  30.  
  31.     fprintf(fpr, "n3  0  0  -2.01\n");
  32.     fprintf(fpr, "*\n");
  33.  
  34.     for (m=0; m<32; m++) 
  35.         fprintf(fpr, "f%03d  n0  n01.m%02d  n01.m%02d\n",i++, m, (m+1)%32);    
  36.  
  37.     for(m=0; m<32; m++) {
  38.         fprintf(fpr, "f%03d  n02.m%02d  n02.m%02d n01.m%02d n01.m%02d\n",i++, m, (m+1)%32, (m+1)%32, m);  
  39.        }
  40.  
  41.     for (m=0; m<32; m++) 
  42.         fprintf(fpr, "f%03d  n3  n02.m%02d  n02.m%02d\n",i++, (m+1)%32, m);    
  43.  
  44.     fprintf(fpr, "*\n");
  45.  
  46.     fclose(fpr);
  47. }
  48.  
  49.  
  50.